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We measure the forward-backward asymmetry in the production of A b and A b baryons as a 
function of rapidity in pp collisions at y/s = 1.96 TeV using 10.4 fb _1 of data collected with the 
DO detector at the Fermilab Tevatron collider. The asymmetry is determined by the preference of 
A b or A*, particles to be produced in the direction of the beam protons or antiprotons, respectively. 
The measured asymmetry integrated over rapidity y in the range 0.1 < |y| < 2.0 is A = 0.04 ± 
0.07 (stat) ± 0.02 (syst). 

PACS numbers: 13.60.Rj, 14.20.Mr 


Hadroproduction of particles carrying a heavy quark 
Q (Q = b,c) proceeds through gluon-gluon fusion or 
quark-antiquark annihilations [l]], followed by hadroniza- 
tion of the heavy quarks. At the parton level of leading- 
order (LO) quantum chromodynamics (QCD), Q and Q 
quarks are produced symmetrically. Next-to-leading or¬ 
der (NLO) QCD effects can introduce a small asymmetry 
of ~ 1% in Q and Q momenta from interfering ampli¬ 
tudes. The liadronization process may also change the 
direction of the particle carrying the quark Q relative to 
the original Q direction and thus generate a significant 
asymmetry. 

There have been few studies of this effect in bottom 
baryon production compared to bottom mesons. Pro¬ 
duction of heavy baryons is sensitive to effects of non- 
perturbative final state interactions of a QCD string con¬ 
necting the b quark and a remnant of the proton. The 
production of the ground-state bottom baryon A b and 
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its antiparticle A b has been recently discussed by Ros- 
ner Q, who proposes the “string drag” mechanism that 
may favor production of A b baryons in the hemisphere 

containing the beam proton, and A b baryons in the an¬ 
tiproton beam hemisphere. In the string drag picture, 
the QCD interaction between a b quark produced in the 
pp collision and the remnant of the proton is described by 
a string with a linear potential. When the string breaks, 
it imparts an impulse to the quark along the beam axis. 
Assuming a string tension of 0.18 GeV 2 , Rosner made an 
approximate prediction for the shift in the particle lon¬ 
gitudinal momentum relative to the axis along the beam 
direction of A p z = 1.4 GeV, resulting in a shift in ra¬ 
pidity of approximately Ay = 1.4 GeV/E, where E is 
the energy of the particle and the rapidity is defined as 
y = ln((.E +p z )/(E — p z ))/ 2. Another possible source of 
asymmetry in A b production is the coalescence of an in¬ 
trinsic b quark at large momentum fraction x in the Fock 
state | uudbb) of the proton with a comoving diquark ud 
from the proton [sj. 

In this article, we present a study of the forward- 
backward production asymmetry of A b and A b baryons 
using the fully reconstructed decay chain A b —> J/if A, 
J/if —> p + p~, A —»• p7T~ , and its charge conjugate. The 
forward (F) category corresponds to a particle (A b or 

A°) sharing valence quark flavors with a beam particle 
with the same sign of rapidity, and the backward ( B) 
category corresponds to the reverse association. In pp 
collisions at DO, we choose the positive z-axis to be in 
the direction of the proton beam, so that the forward di¬ 
rection corresponds to a A b particle emitted with y > 0 
or a A b particle emitted at y < 0. In pp collisions, A b 
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particles are assigned to the forward category and A fc 
particles to the backward category. To facilitate a com¬ 
parison with existing measurements, we present the ra¬ 
tio of the backward to forward production cross sections, 
R = <j(B)/cr(F), and the forward-backward asymmetry, 
A = (a(F) — cr(B))/(a(F) + cr(B)), as functions of the 
rapidity y. The data sample corresponds to an integrated 
luminosity of 10.4 fb _1 collected with the DO detector in 
pp collisions at yfs =1.96 TeV at the Fermilab Tevatron 
collider. 

Using the same data set, the DO experiment has stud¬ 
ied the forward-backward asymmetry in the production 
of B± mesons, observing no rapidity dependence Q. The 
measured forward-backward asymmetry in the produc¬ 
tion of B± mesons, where the forward category corre¬ 
sponds to B~ mesons produced at y > 0 and B + mesons 
produced at y < 0, is Afb(B ± ) = [—0.24 ± 0.41 (stat) ± 
0.19 (syst)]%, integrated over rapidity. 

The DO detector consists of a central tracking system, 
calorimeters, and muon detectors [bij. The central track¬ 
ing system comprises a silicon microstrip tracker and a 
central fiber tracker, both located inside a 1.9 T super¬ 
conducting solenoidal magnet. The tracking system is 
designed to optimize tracking and vertexing for pseudo¬ 
rapidities \p\ < 3, where p = — ln[tan(0/2)], and 9 is 
the polar angle with respect to the proton beam direc¬ 
tion. The tracking system can reconstruct the primary 
pp interaction vertex for interactions with at least three 
secondary tracks with a precision of ~ 35 pm (~ 90 pm) 
in the plane transverse to (along) the beam direction. 
The muon detector, positioned outside the calorimeter, 
consists of a central muon system covering the pseudo¬ 
rapidity region of |7/| < 1 and a forward muon system 
covering the pseudorapidity region of 1 < \p\ < 2. Both 
central and forward systems consist of a layer of drift 
tubes and scintillators inside 1.8 T iron toroidal mag¬ 
nets and two similar layers outside the toroids Q. The 
toroid and solenoid magnet polarities were periodically 
reversed, allowing for a cancellation of first-order effects 
related to a possible instrumental asymmetry. 

Candidate events are required to include a pair of op¬ 
positely charged muons. At least one muon is required to 
be detected in the muon chambers in front of and behind 
a toroid magnet. The other muon may be detected only 
in front of the toroid or as a minimum ionizing particle 
in the calorimeter. Each muon candidate is required to 
match a track found in the central tracking system. 

To form A° and A° candidates, muon pairs in the in¬ 
variant mass range 2.9 < M{p + p~) < 3.3 GeV, con¬ 
sistent with a J/ip meson decay, are combined with A 
baryon candidates. The A candidates are formed from 
pairs of oppositely charged tracks originating from a 
common vertex, consistent with a decay A —> pn~ or 
A — > pir + . The charged particle with the higher momen¬ 
tum is assigned the proton mass. A previous analysis 
has shown that the misassignment of the proton track 
is negligible [7[. The A candidate is required to have 
an invariant mass between 1.107 and 1.125 GeV and a 


transverse momentum greater than 1.8 GeV. The sep¬ 
aration of the A decay vertex from the primary vertex 
in the transverse plane must be between 0.5 and 25 cm. 
A kinematic fit of the parameters of tracks forming the 
A° candidate is performed by constraining the dimuon 
invariant mass to the world-average J/i/j mass Q, and 
constraining the J/ipA system to originate from a com¬ 
mon decay vertex. The modified track parameters are 
used in the calculation of the A^ invariant mass. We 
require 5.0 < M(J/tpA) < 6.2 GeV. 

To suppress the large background from prompt J/ij) 
production, we require a significant separation of the A® 
decay vertex from the primary vertex. To reconstruct the 
primary vertex, tracks are selected that do not belong to 
the A[] decay. We constrain the transverse position of 
the primary vertex to the average beam location in the 
transverse plane. We define the signed decay length of a 
A° baryon, L xy , as the vector pointing from the primary 
vertex to the A[[ decay vertex projected on the direction 
of the A[[ transverse momentum pr- We require L xy to 
be greater than three times its uncertainty. 

The mass distributions for A° candidates in the range 
0.5 < \y\ < 1.0 in the forward and backward cate¬ 
gories are shown in Fig. [l] Binned maximum-likelihood 
fits of a Gaussian signal function and a second-order 
Chebyshev polynomial for the background yield a for¬ 
ward (backward) signal with a mean mass of M(A°) = 
5618.1±4.3 MeV (5619.9±4.7 MeV), consistent with the 
world-average A° mass (sf. The width depends on y and 
varies between about 30 and 50 MeV. The average re¬ 
constructed pr of A° candidates is (pr) = 9.9 GeV after 
background subtraction. 

The production rates of forward and backward A® and 

A° baryons are extracted from fits to the invariant mass 
distributions of forward and backward candidates in four 
rapidity bins in the range 0.1 < \y\ < 2.0, as defined in 
Table [TJ We reject the region \y\ < 0.1 where the asymme¬ 
try may be diluted by forward-backward migration due 
to the finite polar angle resolution [ 4 ] . 

Samples of fully simulated Monte Carlo (MC) signal 
events are obtained at LO with PYTHIA pj and at NLO 
with MC@NLO 0 , using the parton distribution func¬ 
tion sets CTEQ6L1 and CTEQ6M1 [ll|, respectively. 
Pythia generates bb quark pairs via direct 2 —> 2 pro¬ 
cesses ( QiPi.gg —> bb) and decays of gauge bosons, as well 
as through flavor excitation processes like bg —» bg , and 
gluon splittings, g —> bb. The event generator MC@NLO 
is interfaced with Herwig EH for parton showering and 
hadronization. After hadronization, bottom hadron de¬ 
cays are simulated with EvtGen 13]. In the simulation, 

the A° and A & baryons are forced to decay to J/tpA, 
J/4> —> P + P~ , using the phase space (PHSP) and vector 
to lepton-lepton (VLL) models in EvtGen. The detec¬ 
tor response is simulated with GEANT3 0 and multiple 
pp interactions (pileup) are modeled by overlaying hits 
from random bunch crossings in data. A MC sample 
generated with pythia, 30 times the number of signal 
events in the data sample, is used to obtain efficiencies 


5 


TABLE I: Efficiencies e, averaged values of background-subtracted transverse momenta (pt), backward and forward fitted 
yields for the signal N(B) and N(F), forward-backward asymmetries A, and cross-section ratios R in four intervals of rapidity. 
Uncertainties on (j>t), N{B) and N(F) are statistical only. Uncertainties on e arise from the statistical precision of the simulated 
event samples. 


\y\ 

e(%) 

(Pt) (GeV) 

N(B) 

N(F) 

A ± (stat) ± (syst) 

R ± (stat) ± (syst) 

0 . 1 - 0.5 

0.70 ±0.01 

10.2 ±0.1 

125 ± 18 

92 ± 17 

-0.15 ±0.11 ±0.03 

1.36 ±0.32 ±0.06 

0 . 5 - 1.0 

1.01 ±0.01 

10.0 ±0.1 

135 ± 19 

154 ±22 

0.07 ±0.10 ±0.02 

0.88 ±0.18 ±0.04 

1 . 0 - 1.5 

0.97 ±0.01 

9.7 ± 0.1 

123 ± 16 

158 ±23 

0.12 ±0.10 ±0.02 

0.78 ±0.15 ±0.04 

1 . 5 - 2.0 

0.32 ±0.01 

9.8 ±0.2 

22 ±9 

33 ± 10 

0.21 ±0.24 ±0.02 

0.67 ±0.34 ±0.03 



FIG. 1: (color online) Invariant mass distribution of —> 

J/ipA and A;, —» J/ipK candidates in the rapidity range 
0.5 < |j/| < 1.0 in the (a) forward and (b) backward cate¬ 
gories. The fit of a Gaussian signal function with a second- 
order Chebyshev polynomial background function is superim¬ 
posed. The vertical lines define the signal region. 


for reconstructing Aj baryons in each of the four rapidity 
intervals shown in Table [I] The A° efficiencies are sup¬ 
pressed by the large transverse momentum requirement 
on the A candidates and by the low reconstruction effi¬ 
ciency for the long-lived A baryon. 

Most of the systematic uncertainties in the production 
cross-sections of A° and A h baryons arise from uncertain- 


TABLE II: Systematic uncertainties (in %) on the measure¬ 
ment of the backward-to-forward ratio R. 


Source 

Uncertainty (%) 

Signal shape 

2 

Background shape 

2 

Detection efficiency 

3 

Total syst. uncertainty 
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ties in the kinematic acceptance and detection efficiency 
of final-state particles and cancel in the measurements 
of the asymmetry A and ratio R. The remaining un¬ 
certainties are due to the signal and background shapes 
assumed in the mass fits and the different efficiencies of 
forward and backward particle reconstruction. The un¬ 
certainty from the signal shape is estimated by compar¬ 
ing the results of the central fits with the results obtained 
when the width parameters for the forward and backward 
categories are constrained to be equal. The sensitivity 
to the background shape is estimated by increasing the 
lower mass requirement to M(J/ipA) >5.2 GeV, thus ex¬ 
cluding the mass range where feed-down from multi-body 
bottom baryon decays may be present. The estimate of 
the uncertainty on the detection efficiency is based on 
the average deviation from unity of the ratio R of recon¬ 
structed events in four rapidity intervals for a sample of 
MC events generated with no asymmetry. Adding the 
uncertainties in quadrature results in a total systematic 
uncertainty of ±4%. The systematic uncertainties are 
summarized in Table HT1 

The fitted signal yields and the resulting forward- 
backward asymmetry A are presented in Table [T| We ob¬ 
serve that there is a weak correlation between rapidity y 
and the averaged value of background-subtracted trans¬ 
verse momentum ( pt ) of A° candidates. The asymmetry 
integrated over \y\, taking into account the rapidity de¬ 
pendent efficiency e, is A = 0.04±0.07 (stat)±0.02 (syst). 

The forward-backward asymmetry as a function of \y\ 
is shown in Fig. [2] There is a wide range of model pre¬ 
dictions for this asymmetry. The “heavy quark recombi¬ 
nation” model [l5j, as shown in Fig. [2] predicts a modest 
asymmetry, reaching ss 2% near \y\ = 2. While PYTHIA 
predicts no asymmetry, the MC@NLO generator inter- 
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FIG. 2: (color online) Measured forward-backward asymme¬ 
try A versus rapidity \y\ compared to predictions of the heavy 
quark recombination model 15] and a simulated effect of the 
longitudinal momentum shift due to beam drag (see Ref. (2] 
and text). The background asymmetry is obtained from J/ipA 
candidates in the A b mass sidebands (uncertainties are small 
compared to the symbol size). Measurements are placed at 
the centers of the rapidity intervals defined in Tabic [I] 


faced with Herwig predicts a large asymmetry, reaching 
100% close to \y\ = 2. Our results are consistent with no 
asymmetry within the large uncertainties, although they 
show a trend of increasing asymmetry with increasing |y| 
that could be interpreted as the effect of the longitudinal 
momentum imparted to a A b or A b particle by the beam 
remnant. Assuming a shift of A p z = 1.4 GeV in the 
particle longitudinal momentum, as estimated by Ros- 
ner Q, we have simulated the effect by adding 1.4 GeV 
(—1.4 GeV) to the A b (A fc ) baryon p z in the generated 
PYTHIA events. As shown in Fig. (2j our result is in a 
good agreement with this prediction. We find our results 
in disagreement with the large asymmetry predicted by 
MC@NLO+Herwig. 

The results for the backward-to-forward ratio R for 
the same rapidity intervals are given in Table U and 
shown in Fig. [3J where we compare with the results 
for the ratio of cross sections, cr(A b )/cr(A b ), for the six 
rapidity bins reported by the CMS Collaboration 13- 
All results are presented as functions of the “rapid¬ 
ity loss”, defined as the difference between the rapid¬ 
ity of the beam, y(beam) = 7.64 (8.92) at the Teva- 
tron (LHC), and the rapidity y of the A b baryon. The 
DO and CMS results are consistent within large uncer¬ 
tainties. Together, they show a trend of R to fall with 
increasing rapidity and decreasing rapidity loss. The 
DO result for the ratio R integrated over rapidity, tak¬ 
ing into account the rapidity dependent efficiency e, is 
R = 0.92 ± 0.12 (stat) ± 0.04 (syst), to be compared with 
the value of R = 1.02 ± 0.07 (stat) ± 0.09 (syst) reported 
by the CMS Collaboration. 

In order to verify that detector effects on R and A 


FIG. 3: Measured ratio of the backward to forward produc¬ 
tion cross sections versus rapidity loss compared to the A b to 
A b production cross section ratio at CMS taken from Table II 
of Ref. 0- Measurements are placed at the centers of their 
rapidity loss ranges. 


are not significant, the analysis was repeated considering 
candidates with y > 0 (or y < 0) only, and A b (or A b ) 
only. Within statistical uncertainties, all results are con¬ 
sistent with each other and with the measurements listed 
in Table [I] Furthermore, as shown in Fig. [2l we find a 
negligible forward-backward asymmetry in the four in¬ 
tervals of rapidity in a sample of background candidates 
obtained from the A b mass sidebands (region above and 
below the A b signal region defined in Fig. [T| with no L xy 
requirement. 

In summary, we have presented a measurement of the 
forward-backward asymmetry in the production of A b 

and A b baryons as a function of rapidity \y\. Together 
with related results from the LHC, the data show a ten¬ 
dency of forward particles that share valence quarks with 
beam remnants, to be emitted at larger values of rapidity, 
corresponding to smaller rapidity loss, than their back¬ 
ward counterparts. The measured ratio of the backward- 
to-forward production rate at the mean transverse mo¬ 
mentum of (pr) = 9.9 GeV, averaged over rapidity in 
the range 0.1 < \y\ < 2.0, is R = 0.92 ± 0.12 (stat) ± 
0.04 (syst). The measured forward-backward asymmetry 
is A = 0.04 ± 0.07 (stat) ± 0.02 (syst). 
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